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Abstract 

As an outgrowth of our investigation of non-regular spaces within the 
context of quantum gravity and non-commutative geometry, we develop a 
graph Hilbert space framework on arbitrary (infinite) graphs and use it to 
study spectral properties of graph- Laplacians and graph-Dirac-operators. 
We define a spectral triplet sharing most of the properties of what Connes 
calls a spectral triple. With the help of this scheme we derive an explicit 
expression for the Connes- distance function on general directed or undi- 
rected graphs. We derive a series of apriori estimates and calculate it for 
a variety of examples of graphs. As a possibly interesting aside, we show 
that the natural setting of approaching such problems may be the frame- 
work of ( non- Jlinear programming or optimization. We compare our results 
(arrived at within our particular framework) with the results of other au- 
thors and show that the seeming differences depend on the use of different 
graph-geometries and/or Dirac operators. 



PACS numbers: 02.30.Sa, 02.30.Tb, 04.60.-m 



1 Introduction 



In recent years we started a programme to reconstruct continuum physics and/or 
mathematics from an underlying more primordial and basically discrete theory 
living on the Planck-scale (cf. to or sort of a "spin-off" various 



problems of a more mathematical and technical flavor emerged which may have 
an interest of their own. Discrete differential geometric concepts were dealt with 
in 0, the theory of random graphs was a central theme of p|, topics of dimen- 
sion theory and fractal geometry were addressed in lump-spaces and random 
metrics in [p^] . 

If one wants to recover the usual (differential) operators (or more generally, 
the concepts of standard functional analysis), being in use in ordinary contin- 
uum physics and mathematics by some sort of limiting process from their dis- 
crete protoforms, which live, on their part, on a relatively disordered discrete 
background like, say, a network, it is reasonable to analyse in a first step these 
discrete counterparts more closely. This will be one of our themes in the follow- 
ing with particular emphasis on discrete Laplace- and Dirac- operators on general 
graphs. In contrast to |]T6[ we now also include arbitrary directed graphs. The 
main thrust goes however into an analysis of metrical concepts on discrete (non- 
commutative) spaces like general graphs, being induced by graph Dirac- operators 
and the Connes- distance functional. 

We note in passing that functional analysis on graphs is both of interest 
in pure and applied mathematics and also in various fields of (mathematical) 
physics. For one, discrete systems have an increasing interest of their own or serve 
as easier to analyse prototypes of their continuum counterparts. To mention a few 
fields of applications: graph theory in general, analysis on (discrete) manifolds, 
lattice or discretized versions of physical models in statistical mechanics and 
quantum field theory, non- commutative geometry, networks, fractal geometry etc. 
From the widely scattered (mathematical and physical) literature we mention 
(possibly) very few sources we are aware of. Some ot them were of relevance for 
our own motivation, some others we came across only recently (see e.g. to 
1^, and |3^, ||14| or [0). Some more literature like e.g. [0 was pointed 



out to us by Mueller-Hoissen; the possible relevance of references to were 
brought to our notice by some unknown referee. Last, but not least, there is the 



vast field of discretized quantum gravity (see e.g. |^ or [^). All this shows 
that the sort of discrete functional analysis we are dealing with in the following, 
is presently a very active field with a lot of different applications. 

For the convenience of the reader we begin with compiling some concepts 
and tools dealing with graph Hilbert spaces which we then use to investigate the 
spectral properties of graph Laplacians and Dirac operators. In a next step we 
study and test concepts and ideas, which arose in the context of non- commutative 
geometry. As we (and others) showed in preceding papers, networks and graphs 
may (or even should) be understood as examples of non- commutative spaces. 
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A currently interesting topic in this field is the investigation of certain distance 
junctionals on "nasty" or non-standard spaces and their mathematical or physical 
"naturalness" . Graphs carry, on the one hand, a natural metric structure given 
by a distance function d{x, y), with x, y two nodes of the graph (see the following 
sections). This fact was already employed by us in e.g. |Q to develop dimensional 
concepts on graphs. Having Connes' concept of distance in noncommutative 
geometry in mind (cf. chapt. VI of it is natural to try to compute it in 
model systems, which means in our context: arbitrary graphs, and compare it 
with the already existing notion of graph distance mentioned above. (We note in 
passing that the calculation of the Connes distance for general graphs turns out 
to be surprisingly complex and leads to perhaps unexpected connections to fields 
of mathematics like e.g. (non-) linear programming or optimization; see the last 
section) . 

Therefore, as one of many possible applications of our formalism, we construct 
a protoform of what Connes calls a spectral triple, that is, a Hilbert space structure 
, a corresponding representation of a certain (function) algebra and a (in our 
framework) natural candidate for a so-called Dirac operator (not to be confused 
with the ordinary Dirac operator of the Dirac equation), which encodes certain 
properties of the graph geometry. This will be done in section ^. 

In the last (and central) section, which deals with the distance concept deriv- 
ing from this spectral triplet we will investigate this concept more closely as far 
as graphs and similar spaces are concerned. In this connection some recent work 
should be mentioned, in which Connes' distance function was analyzed in certain 
simple models like e.g. one-dimensional lattices (|^-[^). These papers already 



show that it is a touchy business to isolate "the" appropriate Dirac operator (after 
all, different Dirac operators are expected to lead to different geometries!) and 
that it is perhaps worthwhile to scrutinize the whole topic in a more systematic 
way. We show in particular that one may choose different Dirac-operators on 
graphs (or rather, different types of graphs over the same node set) which may 
lead to different results for e.g. the corresponding Connes- distance. 

The problem of finding suitable metrics on "non-standard" spaces is a par- 
ticularly interesting research topic of its own, presently pursued by quite a few 
people (see the papers by Rieffel or Weaver [Q,|jl3| and the references mentioned 
therein). Another, earlier (and important) source is [l^. We recently extended 



the investigation of metric structures to socalled lump spaces and probabilistic 
metric spaces (see and |3^) and employed it in the general context of quan- 
tum gravity (cf. also 0). 

Remark: When we wrote an earlier draft of this paper we were unaware of the 



content of |T^. It happened only recently that we realized that various of the 
results we derived in connection with Dirac-operators and the Connes-distance 
can already be found in and we try to take care of this fact in the follow- 
ing. Both the technical approach and the motivation are, however, not entirely 
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identical. The interested reader may consult the electronic version of our earlier 
draft (PI). 



The same applies to paper 10 1 where some of the Hilbert space methods 



were developed we later rederived in |16| being unaware of the prior results. As 



a consequence we drop most of the technical steps leading to that part of our 
previous results overlapping with corresponding parts in the above mentioned 
papers and refer, for simplicity to |jl6| where the interested reader can find more 
details. 

2 A brief Survey of Differential and Operator 
Calculus on Graphs and Graph 
Hilbert Spaces 

We give a brief survey of certain concepts and tools needed in the following 
analysis. While our framework may deviate at various places from the more 



traditional one, employed in e.g. algebraic graph theory (see e.g. 0,0 or pTI), 
this is mainly done for reasons of greater mathematical flexibility and generality 
and, on the other side, possible physical applications (a case in point being the 
analysis of non- commutative spaces). Most of the technical tools which are not 
defined in detail in the following have been introduced in section 3 of 

2.1 Simple (Symmetric) Graphs 

For simplicity we assume the graph to be connected and locally finite (no el- 
ementary loops and multiedges, whereas these could easily be incorporated in 
the framework), i.e. each node (or vertex) is incident with only a finite num- 
ber of edges (or bonds). To avoid operator domain problems we usually make 
the even stronger assumption that the vertex- or node degree, Vi, {i labelling the 
nodes), is globally bounded but this restriction is frequently not really necessary. 
Furthermore, it has turned out to be algebraically advantageous to identify an 
(undirected) labelled graph with a directed graph having two oppositely directed 
edges for each undirected edge, the directed edge, pointing from node Ui to node 
Uk, being denoted by dik, the oppositely directed edge by dki and the undirected 
(but orientable) edge by their superposition 

bik = dik — dki = —bki (1) 
{dik and dki are treated as independent basis vectors; cf. ^ or HTBI). 



As the elementary building blocks of our graph Hilbert spaces we take {nj} 
and {dik} as basis elements of a certain hierarchy of Hilbert spaces over, say, C 
with scalar product induced by 

{rii\nk) = 6ik {dik\dim) = Su ■ 5km (2) 
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This definition implies {bik\bik) = 2. 

Definition 2.1 (Vertex-, Edge Hilbert Space) The Hilbert spaces Hq, 
(a for antisymmetric) and Hi consist of the formal sums 

f-=^fini g := y^jikdjk with gik = -gki and g' ■.= ^gikdik (3) 



fi,gik ranging over a certain given field like e.g. C (sometimes only rings like 
e.g. Z are admitted; then we are dealing only with modules). We evidently have 
HI C Hi. 

Remark: One could continue this row of vector spaces in ways which are common 
practice in, say, algebraic topology ( see |l|] sections 3.1 and 3.2). In this context 
they are frequently called chain complexes (see also |^). On the other hand, the 



above vector spaces could as well be viewed as discrete function spaces over the 
node-, bond sei with ni,dik now representing the elementary indicator functions. 



Proceding in this spirit we can now introduce two linear maps between Hq, Hi 
extending the usual boundary- and coboundary map. On the basis elements they 
act as follows: 

6 : dik^ Uk hence bik Uk - Ui (5) 
d-rii^ ^(dfei - dik) = ^ bki (6) 

k k 

and linearly extended. That is, 6 maps the directed bonds dik onto the terminal 
node and b^k onto its (oriented) boundary, while d maps the node Ui onto the sum 
of the ingoing directed bonds minus the sum of the outgoing directed bonds or 
on the sum of oriented ingoing bonds bki- 

As was shown in [|l| (we later realized that the same definition was already 
employed in []T3|), these definitions lead in fact to a kind of discrete differential 
calculus on Hq,Hi, that is we have 

df = diJ2 f^^d = Y.^fk - fddik (7) 

k,i 

Combining now the operators 5 and rf, we can construct the canonical graph 
Laplacian. On the vertex space it reads: 

= - EE - ■ = - - /^))^^ =■ -^f (8) 

i k i k 
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where Vi denotes the node degree or valency defined above and the fc-sum extends 
over the nodes adjacent to n^. 



Remark: Note that there exist several variants in the hterature (see e.g. or 
0). Furthermore, many mathematicians employ a diff'erent sign- convent ion. We 
stick in the following to the convention being in use in the mathematical-physics 
literature where —A is the positive(!) operator. 

This graph Laplacian is intimately connected with another important object, 
employed in algebraic graph theory, i.e. the adjacency matrix, A, of a graph, 
its entries, Uik, having the value one if the nodes ni,nk are connected by a bond 
and are zero elsewhere. If the graph is undirected (but orientable), the relation 
between ni,nk is symmetric, i.e. 

Oifc = 1 =^ aki = 1 (9) 

This has the obvious consequence that in case the graph is simple and undirected, 
y4 is a symmetric matrix with zero diagonal elements. 

Remark: More general A's occur if more general graphs are admitted (e.g. general 
multigraphs). 

With our definition of A it holds: 

A = A-V (10) 

where V is the diagonal degree matrix, having Vi as diagonal entries. (Note that 
the other sign-convention would lead to A = V — A). 

To avoid domain problems we assume from now on that the node degree, Vi, 
is uniformly hounded on the graph G, i.e. 

Vi < Vmax < OO (11) 

Defining di^2 as 

di,2 ■ rii-^^dki, ^dik (12) 
respectively and linearly extended, we get 

d = di — d2 (13) 
Similarly we make the identification 6 =: 6i with 

^1,2 : dik Uk, Ui (14) 

It is noteworthy (but actually not surprising) that Vi < Vmax implies that 
all the above operators are hounded (in contrast to their continuous counter- 
parts, which are typically unbounded). Taking this for granted at the moment, 
a straightforward analysis yields the following relations: 
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Lemma 2.2 

1. The adjoint d* of d with respect to the spaces Hq,Hi is 26 

2. On the other side we have for the natural extension of d, 6 to the larger 
space Hi: 

Si = {diY , 52 = {d2T (15) 

hence 

{Si - S2) = (di - d2y = d*^2S = 2Si (16) 

3. Furthermore it holds 

dl ■ di = 61 ■ di = dl ■ d2 = V : rii Vitii (17) 
(with V the vertex degree matrix) 

dl ■ d2 = Si ■ d2 = dl ■ di = S2 ■ di = A : Ui -* ^ Uk (18) 

Similar geometric properties of the graph are encoded in the products coming 
in reversed order. 

(Here and in the following k — i means summation over the first index and runs 
through the set of labels of nodes directly connected with n^). 

That and how d, d* encode certain geometric information about the graph 
can be seen from the following domain- and range-properties (cf. fl^, for corre- 
sponding results in the more traditional approach see [p,p.24ff). 

Theorem 2.3 Let the graph be connected and finite, \V\ = n, then 

dim{Rg{d*)) =n-l (19) 

dim{Ker{d*)) = ^Vi- (n-l) (20) 

i 

With dim{Hi) = dim{H^) = 1/2 ■ dim{Hi) we have 

codim{Ker{d*)) = dim{Rg{d)) =n-l (21) 
We see that both Rg{d*) and Rg{d) have the same dimension {n — 1). 
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Remark 2.4 In case the graph has, say, c components, the above results are 
altered in an obvious way; we have for example 

dim{Rg{d*)) = n-c (22) 

In the literature Ker{d*) is called (for obvious reasons) the cycle subspace (cf 
e.g. [§). On the antisymmetric subspace we have d* = 26 and 6{bik) = 
rik — rii. Choosing now a cycle, given by its sequence of consecutive vertices 
nj,, . . . , ni^^Ui^^^ := n^,, we have 

that is, vectors of this kind lie in the kernel of d* 

We will now provide quantitative lower and upper bounds for the respective 
norms of the occurring operators. We have: 

iM/ir = E K^*^ - = ■ i^^i' + E^'^ ■ i^^i' - E(^^ + 7^fk) 

ik i k ij^k 

= 2-E^^I/^r-2-E^/^ (24) 

i i^k 

which can be written as: 

IWr = m\Vf) - {f\Af)) = (/I - 2A/) (25) 

and shows the close relationship of the norm of d with the expectation values of 
the adjacency and degree matrix or the graph Laplacian. That is, norm estimates 
for, say, d, derive in a natural manner from the corresponding estimates for A or 
-A. With 

\\dfr = if\d*df) = {f\-2Af) (26) 

we have 

0<d*d = -2A and \\df = sup (/| - 2A/) = || - 2A|| (27) 

11/11=1 

Furthermore via 

< sup (/I - 2A/) < 2v^a. + 2 sup I < f\Af > \ (28) 

11/11=1 ll/il=i 

we get 

II - All < Vmax + \\A\\ (29) 
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Remark: We want to mention that we are using the usual operator norm also 
for matrices (in contrast to most of the matrix literature), which is also called 
the spectral norm. It is unique in so far as it coincides with the so-called spectral 



radius (cf. e.g. p9[ or |3y]), that is 

\\A\\ := sup{|A|; A G spectr{A)} (30) 

In a first step we give upper and lower bounds for the operator norm of the 
adjacency matrix, A, both in the finite- and infinite-dimensional case. There 
are various proofs available of a varying degree of generality (see e.g. [10 , [1? 



or ||T^) to which we refer the reader. In the following we give only the final 
results. Note however that the transition from finite to infinite graphs is far from 
straightforward as in some of the necessary technical steps entirely new methods 
are needed. 

Theorem 2.5 (Norm of A) With the adjacency matrix A finite or infinite and 
a finite Vmax we have the following result (a certain fixed labelling of the nodes 
being assumed): 

n 

limsup n^^ ■ ''^^vf^^ < limsup \\An\\ = \\A\\ = sup{|A|; A G spectr{A)} < Vmax 

(31) 



i=l 



Here A^ are the adjacency matrices for the induced subgraphs, living over the 
first n labelled nodes, f -"^ is the corresponding induced (and n-dependent) node 
degree. 

As a byproduct we have the important lemma 

Lemma 2.6 The adjacency matrices, An, converge strongly to A and we have 
in particular \\An\\ / 



(For a proof of the latter result see [T^ or |16| 



Remark 2.7 To prove strong convergence of operators is of some relevance for 
the limit behavior of spectral properties of the operators An, A. That is (cf. 
e.g. section VIII. 7), we have in that case (An, A self adjoint and uniformly 
bounded) An ^ A in strong resolvent sense, which implies that the spectrum of 
the limit operator. A, cannot suddenly expand, i.e. 

A G spec{A) ^ 3 A„ G spec{An) with A„ — A (32) 

and for a,b ^ speCpp{A) 

Pia,b){An) P{a,b)iA) strongly (33) 
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To test the effectiveness of tlie upper and lower bounds given above, we apply 



them to a non-trivial model recently discussed in i.e. the infinite binary tree 
with root no where vq is two and Vi equals three for i ^ 0. The authors show 
(among other things) that the spectrum consists of the interval [— 2v^, 2-\/2], 
i.e. ||y4|| = 2a/2. Vmax is three, we have to calculate limsupl/n ■ -^^^ 
simplicity we choose a subsequence so that n := n{N) with denoting the A^-th 
level (consisting of 2^ nodes) of the tree starting from the root no. Note that in 
the corresponding induced subgraph the boundary nodes sitting in the A^-th 
level have only node degree one with respect to Gat but three viewed as nodes in 
the full tree. 
We then have 



n 



Hence 



N n{N) N-1 N 

^ 2\ ^ t;, = 2 + 3 • ^ 2^ + 2^ = 3 ■ J] 2*= - 2 ■ 2^ - 1 (34) 

k=l i=0 k=l k=0 



n(N) N 



lim l/n(N) V = 3 - 2 lim( V 2'"^^)-^ = 2 (35) 

n(N) ' ' /-^ N 

^ ' i=0 

That is, our general estimate imply 2 < ||y4|| < 3, which is not so bad. 



2.2 Arbitrary Directed Graphs 

If we deal with general directed graphs we have both ingoing and outgoing edges 
at each node but in general they do no longer occur in a symmetric way. But nev- 
ertheless, most of our concepts and tools, developed in the foregoing subsection, 
do still exist. The definitions of d and d* are unaltered. As each edge, djfc, is an 
outgoing edge for node Ui and an ingoing edge for node Uk, the same expression 
holds for df, i.e. 

df = ^ifi - fk)dik with / = ^ fiUi (36) 

i,k i 

(the sum of course only extends over those directed edges which do exist in the 
directed graph; in particular each edge, dik, is only counted once in the sum as an 
outgoing edge with respect to the label i.) Furthermore the notions of di^2 and 
^1^2 remain the same, mapping nodes or edges on ingoing edges, outgoing nodes 
and vice versa. We again have 

^1,2 = c^i 2 (5i — ^2) = d* = {di — ^2) (37) 

We can now calculate dl ■ di and ■ ^2 and get: 

dl ■ d,in,) = vr ■ m , d; ■ d2ini) = ■ (38) 
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with the in-, out- degree of the node n^. In the same way we can calculate 

5i ■ d2{ni) = ^ rifc 

k—i,out 

This yields 

d*d = {6i — 62) ■ {di — c/2) 

where the occurring operators on the rhs are the in-, out-vertex degree matrices, 
in-, out- adjacency matrices respectively. 

In general the individual in-, out-adjacency matrices are non-symmetric. Their 
sum, however, is symmetric (and is, in the case of a symmetric graph, twice the 
adjacency matrix of the undirected graph, i.e. 2A). One can now again define a 
(positive, i.e. symmetric) Laplace operator for a non-symmetric (directed) graph, 
that is: 

Conclusion 2.8 

-A := d*d = {V'" + 1/""*) - (A*" + A°"*) =: Vd - Ad (41) 
(Note the now missing factor two in front of A!) 

3 The Spectral Triplet on a general (directed or 
undirected) Graph 

We begin by making some remarks on various concepts, being in use in the 
more recent literature. We note that our version of a Dirac operator (defined 
below) intertwines node- vectors and bond- vectors while in other examples it maps 
node- to node-functions. Our bond-functions have (in some sense) the character 
of cotangential-vectors, while in other approaches derivatives of functions are 
interpreted as tangent-vectors. In our view, the latter formalism is only effective 
in certain classes of highly regular models (like e.g. lattices) where one has kind of 
global directions and will become cumbersome for general graphs. We developed 
this latter approach a little bit in section 3.3 of ||l| and showed how these cotangent 
and tangent vectors can be mapped into each other. We think that, on the 
other side, our framework is more flexible in the general case. This holds in 
particular for our Dirac operator, which encodes certain geometric properties of 
the underlying discrete "manifold" . 

The Hilbert space we will use in the following is 

H = Ho®Hi (42) 



k—i,in 



(39) 



{V' + \/°"*) - {A"' + A°"*) (40) 
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The natural representation of the function algebra (consisting of the bounded 
node functions) 

{f-feC^i.e. sup|/,|<oo} (43) 

i 

on H by bounded operators is given by: 

Ho-- /■/' = ^/./;-n,for/'Gifo (44) 



^1 : / ■ 5Z fidikdik (45) 

From previous work ([|l|) we know that Hi carries also a right-module structure, 
given by: 

gikdik ■ / •= X] ■ ^^'^ ^'^^^ 

(For convenience we do not distinguish notationally between elements of and 
their Hilbert space representations). 

Remark: The same bi-module structure and the Dirac operator defined below 
were already employed in ||T^, p.414ff. 

An important object in various areas of modern analysis on manifolds or in 
Connes' approach to noncommutative geometry is the so-called Dirac operator D 
(or rather, a certain version or variant of its classical counterpart; for the wider 
context see e.g. 0] or |Q to [^). As D we will take in our context the operator: 



acting on 



with 

= (<5i - ^2) (49) 

Note however, that there may exist in general several possibilities to choose such 
an operator. On the other hand, we consider our personal choice to be very 
natural from a geometrical point of view. 
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Lemma 3.1 There exists in our scheme a natural chirality- or grading operator, 
X and an antilinear involution, J. given by 



1 
-1 



with 



and 



[X,^] = x-D + D-x = 



J : 



(50) 



(51) 



(52) 



so that 



J-f-J = f (53) 
These are some of the ingredients which estabhsh what Connes calls a spectral 



triple (cf. e.g. [23] or [24]). We do not want, however, to introduce the full 
machinery at the moment as our scheme has an independent geometric meaning 
of its own. Note in particular what we are saying below about (non) compactness 
of various operators in observation |3.3| . 

Definition 3.2 (Spectral Triplets) As spectral triplet on a general graph we 
take 



(54) 



In our general framework we got in a relatively straightforward manner a 
Hilbert space being the direct sum of the node space (a function space) and the 
bond space (resembling the set of cotangent vectors) and a Dirac operator which 
emerged naturally as kind of a square root of the Laplacian. 

On the other side, if one studies simple models as e.g. in [^ to |2^, other 
choices are possible. In [^,[^, where the one-dimensional lattice was studied, 
the symmetric difference operator was taken as Dirac operator. In [^ the one- 
dimensional lattice was assumed to be directed (i.e., in our notation, only (ij^j+i 
were present) and the Dirac operator was defined (somewhat adhoc) as a certain 
self adjoint "doubling" of the (one-sided, i.e. non-symmetric) adjacency matrix. 
As we will show below, this latter model fits naturally in our general approach 
which includes both directed and undirected graphs. All these Dirac operators 
are different and it is hence no wonder that they lead to different consequences 
(see below). It is our opinion that, in the end, an appropriate choice has to be 
dictated by physical intuition. Nevertheless, this apparent non-uniqueness should 
be studied more carefully. 
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As can be seen from the above, the connection with the graph Laplacian is 
relatively close since for e.g. a symmetric graph we have: 

^2 _ f d*d \ 
\ dd* J 

and 

d*d = -2A (56) 

dd* is the corresponding object on Hi. (In the vector analysis of the continuum 
the two entries correspond to divgrad, graddiv respectively ). 

In the original approach of Connes great emphasis was laid on the compactness 
of operators like the inverse of the Dirac operator and it is part of the general 
definition of a spectral triple. As discrete spaces of the kind we are studying 
are non-trivial examples of non-commutative spaces, it is interesting that we can 
easily test whether this assumption is fulfilled in our particular setting. As may 
be expected, for graphs with globally bounded node degree, we have the following 
result: 

Observation 3.3 Note that all our operators are bounded, the Hilbert space is 
(in general) infinite dimensional, hence there is no chance to have e.g. {D — z)~^ 
or [D"^ — z)~^ compact. At the moment we are sceptical whether this latter phe- 
nomenon dissappears generically if the vertex degree is allowed to become infinite. 
There are some results on spectra of random graphs which seem to have a certain 
bearing on this problem (cf. e.g. P^/j. 

In the next sections we introduce and calculate the socalled Connes-distance 
functional and compare it, among other things, with the ordinary graph distance. 
In doing this we have to calculate the commutator [D, /] applied to an element 
f'eHo. We have: 

{d■f)f = J^{fkfL-md^k (57) 

ik 

if-d)f' = Y.f\{f',-f:)d,k (58) 

ik 

hence 

[DJ]r = J2(f>'-f^)fkd,, (59) 

ik 

On the other side the right-module structure allows us to define df as an operator 
on Hq via: 

df-f= {Y,{fk - fdd^k) ■ {J2 f'k^k) = Y.(f>^ - f^)fkd.'^ = [D, f]f' (60) 

ik k ik 
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(55) 



In a next step we define df as operator on Hi which is not as natural as on Hq. 
We define: 

dflHi ■ dik ifi - fkW (61) 
and hnearly extended. A short calculation shows 

dfU = -{df\H,y = [d*J] (62) 

with 

[d%f]g = d*{f-g)-fd*g {g e H{) (63) 
This then has the following desirable consequence: 

Conclusion 3.4 With the above definitions the representation of df on H is 
given by 



dfUA ^ ( -{df\H,y 

1^ \df\H, ; 



(64) 



df\H=[u\^ ^^^y^HtA/] (65) 



and it immediately follows 

4 The Connes-Distance Function on Graphs 

From the general theory of operators on Hilbert spaces we know that: 

||r|| = ||T*|| (66) 

Hence 

Lemma 4.1 

IIK/]II = IIK/]II = IIM*,/]|| (67) 

and 

ll[A/]|| = |IK/]|| (68) 



Proof: The left part of (|67D is shown below and is a consequence of formula (|7^ ; 



the right identity follows from (pBI). With 



^ - ( ^ ) (69) 
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and Ti := [d, /], T2 := [d* , /], the norm of [D, /] is: 

\\[D, /]f = sup{||Tixf + \\T,yf- ||xf + ||i/f = 1} (70) 

Normalizing now x,y to = = 1 and representing a general normalized 
vector X as: 

X = Ax + /xy , A, /i > and + /i^ = 1 (71) 

we get: 

II [D, /] f = supjA^llTixf + ^?\\T,yr■, \\x\\ = \\y\\ = + = 1} (72) 

where now x, 7/ can be varied independently of A, /i in their respective admissible 
sets, hence: 

||[D,/]f = sup{AlTif + /i2||T2f} = ||Tif □ (73) 



(as a consequence of equation 

We see that in calculating ||[Z},/]|| we can restrict ourselves to the simpler 
expression ||[(i, /]||. We infer from the above calculations {x G Hq): 

i j = l 

and the corresponding expression for a directed graph with Vi replaced by f 
Abbreviating 

J2\fk,-f^\'=■■^^>0 (75) 

i=i 

and calling the supremum over i a^, it follows: 

\\df ■ xjp = as ■ (^tti/tts ■ Ixip) < as (76) 

i 

for \\x\\^ = J2i l^jP = 1- 

On the other side, choosing an appropriate sequence of normalized basis vec- 
tors so that the corresponding a^, converge to we get: 

IM/-e.f^a, (77) 

We hence have 
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Theorem 4.2 



||[A/]||=sup(^|/,^,-/,P)V2, sup(5^|/,^.-/,p)V2 ^,,p,,t,^,ly (78) 

' J=l ' j=l 

The Cannes- distance functional between two nodes, n, n', is now defined as 
foUows: 

Definition 4.3 (Connes-distance function) 

dzstcin,n') := sup{|/„, - \\[D,f]\\ = \\df\\ < 1} (79) 



We would like to note that Davies in [|15[ introduced several metrics on graphs, 
which have been motivated, as he remarks, by his study of heat kernels on Rie- 
mannian manifolds. What he calls metric d^, is related to the rhs of the equation 
in theorem ^]2| (he uses slightly different Hilbert spaces). He then shows in a 
longer proof that this metric is identical to another one, d^, which corresponds 
to the Ihs in the above theorem. In our approach, on the other side, the content 



of theorem |4.2| is derived in a relatively transparent and straightforward way. 



Remark 4.4 It is easy to prove that this defines a metric on the graph. 

Corollary 4.5 It is sufficient to vary only over the set {/; \\df\\ = 1}. 
Proof: This follows from 

\fk-fi\=c-\fk/c-fi/c\; c=\\df\\ (80) 

and 

|M(//c)l|=c-^lM/||=l (81) 
with c < 1 in our case. □ 

In general it turns out to be a nontrivial task to calculate this distance on 
an arbitrary graph as the nature of the above constraint is quite subtle . The 
underlying reason is that the constraint is, in some sense, inherently non-local. 
As / is a function, the difference, fn' — fn, has to be the same independently of the 
path we follow, connecting n' and n. On the other side, in a typical optimization 
process one usually deals with the individual jumps, fk — fi, between neighboring 
points along some path. It is then not at all clear that these special choices of 
jumps along such a path can be extended to a global function without violating 



the overall constraint on the expression in theorem |4.2| . Nevertheless we think the 
above closed form is a solid starting point for the calculation of distc on various 
classes of graphs or lattices. We illustrate this by proving some apriori estimates 
concerning this distance function and by evaluating it for some examples. 
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4.1 Some General Estimates 

Having an admissible function / so that supj(X]fc=i l/fc — /iP)^^^ < 1, this imphes 
that, taking a minimal path 7 from, say, n to n', the jumps \f,y+i — f„\ between 
neighboring nodes along the path have to fulfill: 

l/.+i-M<l (82) 

and, a fortiori, have to be strictly smaller than 1 in the general situation. 

On the other side the Connes distance can only become identical to the ordi- 
nary distance d{n, n') if there exist a sequence of admissible node functions with 
all these jumps approaching the value 1 along such a path, which is however 
impossible in general as can be seen from the structure of the constraint on the 
expression in theorem ^]2| . Only in this case one gets: 

I Y.^Ui - /.) I 5^ 1 = length{^) (83) 
7 7 

We formulate this observation as follows: 

Lemma 4.6 (Connes-distance) Within our general scheme one has the fol- 
lowing inequality 

distc{n,n') < d{n,n') (84) 

By the same token one can prove that distc between two nodes is hounded by 
the corresponding Connes-distance calculated for the (one-dimensional) sub-graph 
formed by a minimal path connecting these nodes, i.e. 

distc{n,n') < distc {rnin. path) {n,n') (85) 

The reason is that one has more admissible functions at ones disposal for a 
subgraph. With G' a connected subgraph of G, the set of admissible function, 
Sg', on G' contains the restrictions of the functions of the corresponding set, Sq, 
belonging to G, as each restriction to G' of a member belonging to Sg lies in 
5*^'. Hence the supremum is in general larger on 5*^'. The distance along such a 
path, on the other side, can be rigorously calculated (see the discussion of some 
examples below) and is for non-neighboring nodes markedly smaller than the 
ordinary graph distance. From what we have said we can also infer the following 
corollary. 

Corollary 4.7 With G' a connected subgraph of G it holds (with n,n' & V G V ) 

distc{n, n'; G) < distc{n, n'; G') (86) 

One can also give sufficient criteria for distc{n,n') < d{n,n'). The cases of 
undirected, directed graphs, respectively, have to be treated a little bit differently. 
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Lemma 4.8 LetG be an undirected graph and^ a minimal path of length, I > 1, 

connecting n, n' . There is at least one node, n* , belonging to 7, having node degree 
>2 (as there are at least two consecutive edges, belonging to j). If 

distc{n, n') — I — d{n, n') (87) 

all the individual jumps along 7 have to be one. But then the corresponding 
function cannot be admissible at n* . Hence 

distc{n,n') < d{n,n') (88) 

Let G now be a directed graph and let there exist two different paths, 7, 7', of 
equal length, I > 1, connecting n,n'. Again there exists a node, n* , on 7 so that 
it is incident with two edges, the one belonging to 7, the other to 7'. Along both 
edges the jumps have to be one and again the admissibility of the corresponding 
function is violated. We again can conclude 

distcin, n') < d{n, n') (89) 

Remark:The latter situation will be discussed below in the example of the directed 
Z^-lattice. 

We remarked above that the calculation of the Connes distance on graphs 
is to a large part a continuation problem for admissible functions, defined on 
subgraphs. Then the following question poses itself. For what classes of graphs 

and/or subgraphs do wc have an equality in the above corollary? We start from 
a given graph, Gq = (Vq, Eq), and then add new nodes and bonds, yielding a new 
graph, G' = (y',E'). We consider two fixed nodes, UojUq in Go- 
Assumption 4.9 We assume that the above process does not create new paths 
between nodes belonging to Gq. In other words, the paths, connecting no and 
are contained in Gq. 

Lemma 4.10 Under this assumption each admissible function on Gq can be ex- 
tended to an admissible function on G' . 

Proof: In a first step we construct the set of nearest neighbors., Vi\Vo in 
relative to Vq. Each new node in Vi\Vo has a unique nearest neighbor in Vq since 
otherwise there would exist a new path between these two nodes lying in Gq. 
With n e 14 \ we extend an admissible function on Gq as follows: 

fn '■= /no ^0 the unique nearest neighbor in Vq (90) 

This extended function is an admissible function on Gi := (Vi, £"1). Note however 
that, by assumption, there do not exist bonds in Ei, connecting nodes in Vi\Vq. 
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We can now continue this process until we arrive at the graph G' . □ 
By the same token we see that 

distc{no, Uq; Gi) = distc{no, Uq] Go) (91) 
This holds at every intermediate step and wc get: 
Lemma 4.11 Under the above assumption we have 

distc{no, Uq-, Gq) = distc{no, Uq; G') (92) 

Corolleiry 4.12 If G is a tree, it holds 

distc{n; n') — distc{n; n'; minimal path) (93) 

Proof: In a tree there exists, by definition, at most one path, connecting two 
nodes. We can take this path as connected subgraph. Go, and make the above 
extension, since G and Gq fulfill the assumption. □ 

The graphs, so constructed are however rather special, consisting, so to speak, of 
a start graph plus some added hair. 

Above we have given sufiicient conditions for 

distc{n, n'; Gq) = distc{n, n'', G) (94) 

with G an extension of Gq and n,n' G Vq. We show now that the emergence of 
too short new paths is representing the obstruction for such a result to hold in 
general. 

So let again Go be a graph and assume the existence of two nodes, n, n', in 
Vq with 

distcin,n';Go) > I eN (95) 

We extend Gq to some G by adding new nodes and edges. We know that for 
admissible functions the elementary jumps along an edge have to fulfill \fi — fk\ < 
1. If there exists a new path, 7, in G, connecting n,n', with 

length{-f) < I (96) 

we can conclude that for each admissible function on G it must hold 

\f{n)-f{n')\<l (97) 

We hence have 
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Lemma 4.13 If two nodes in Go have distc{n,n';GQ) > I and if there exists a 
path, 7; in G, connecting n, n' and having length I <l, it necessarily holds that 

distc{n, n'; G) < I < distc{n, n'; Gq) (98) 

Up to now we have derived upper bounds on distc relative to distc on sub- 
graphs or the canonical graph distance d{n,n'). In the following we will derive 
a quite efficient lower bound. This is done by defining a particular admissible 
function, depending on an arbitrary base node, uq. That is, we fix an arbitrary 
node, dubbed tiq, and take as admissible function the canonical distance, divided 
by the local vertex degree: 

/no(n) := M-'/^ ■ d{no,n) , A„(n) := {v^^'Y'/^ ■ d{no,n) , A„(no) = (99) 

for undirected, directed graphs, respectively. 
From our general results we have 

\\df\\ ^ sup \^V,^-f,\'j (100) 

or Vi replaced by vf'^^K Inserting the above particular function we get 

\\df\\ < 1 (101) 

as each term, \fk^ — fi\, is either zero or one (depending of whether the distance 
to the base point remains constant or changes by ±1). 

Lemma 4.14 The functions, fnoin), rig an arbitrary node in G, are admissible. 
With n, n' two arbitrary nodes in G we take n as base point, Uq, and have 

fnoin') - fnoin) = fnoin') = (vn')-'^' ' d{n' , n) (102) 

(as fno{n) = /no(^o) = 0). As distc is the supremum over admissible functions, 
we get: 

Theorem 4.15 

distcin, n') > {v^n,n'))~^^'' ■ d{n, n') (103) 
with V(n,n') the minimum of the (out-) vertex degrees at n,n' respectively. 

Note that one can of course either choose n or n' as base point in the definition 
of the above admissible function. 
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4.2 Examples 

The general results derived above should be compared with the results in e.g. 
2^ to 1^3. Choosing the symmetric difference operator as "Dirac operator" 



in the case of the one- dimensional lattice the authors in |25, Eg] got a distance 



which is strictly greater than the ordinary distance but their choice does not 



fulfill the above natural constraint given in Theorem ^4.2| . Note in particular that 
our operator ci is a map from node- to bond-functions which is not the case in 
these examples. In ||2^ the authors employed a symmetric doubling of the non- 
symmetric adjacency matrix of the one-dimensional directed lattice, as Dirac 
operator. With v^"^^^ = 1 in this example, our above general estimate yields 

d{n,n') < distc{n,n') < d{n,n') =^ distc{n,n') = d{n,n) (104) 

that is, we get the same result for our Dirac operator as for the choice made in 



We want to close this paper with the discussion of several examples which show 
that, in general, it is quite a non-trivial task to calculate distc- The first one is 
a simple warm-up exercise, the second one is the one-dimensional non-directed 
lattice, 7} , discussed also by some of the authors mentioned above (treated how- 
ever within their own schemes) and is not so simple. The last one is the directed 
Z^-lattice, which we do not solve in closed form, but we provide several estimates. 

The technique used in approaching some of the problems may be interesting in 
general. It turns out that the proper mathematical context, to which our strategy 
does belong, is the field of (non-) linear programming or optimization (see e.g. [BB 



or any other related textbook). This can be inferred from the structure of the 
constraints we get. This means that the techniques developed in this field may 
perhaps be of use in solving such intricate problems. 

Example 1: The square with vertices and edges: 

Xi — X2 — — Xi — Xi (105) 

Let us calculate the Connes-distance between xi and 0:3. As the sup is taken 
over functions, the summation over elementary jumps is (or rather: has to be) 
path- independent (this represents a subtle constraint for practical calculations). 
It is an easy exercise to see that the sup can be found in the class where the two 
paths between xi,X3 have the valuations (1 > a > 0): 

Xi — X2 : a , X2 — X3 : (1 — a^)^^^ (106) 



Xi 



X4 : (1 — a'^y^'^ , X4 — X3 : a (107) 



Hence one has to find supQ<„<]^(a + Vl — a^). Setting the derivative with respect 
to a to zero one gets a = \/l/2. That is: 
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Example 4.16 (Connes-distance on a square) 

distc{xi, X3) = \pl < 2 = d{x\, X3) 



(108) 



Remark 4.17 As Vi = 2, our apriori estimate in theorem [^.i^ is saturated as 

distcixi,X2) > (2)-^/2 ■2 = V2 (109) 
The next example is considerably more complicated. 
Example 2: The undirected one- dimensional lattice: 

The nodes are numbered by Z. We want to calculate distc{0,n) within our 
general framework. The calculation will be done in two main steps. In the first 
part we make the (in principle quite complicated) optimization process more 
accessible. For the sake of brevity we state without proof that it is sufficient to 
discuss real monotonely increasing functions with 

I / [n) tor k > n 

and we write 

k 

f[k) = f{0) + '^h, for 0<k<nhi>0 (111) 

i=l 

The above optimization process then reads: 

Observation 4.18 Find sup ^2^=1 ^« under the constraint 

hl<l,hl + hl<l,... ,hl + hl_, <l,hl<l (112) 

The simplifying idea is now the following. Let h := (/ii)"=i be an admissible 
sequence with all hf_^_^ + hf < 1. We can then find another admissible sequence 
h' with 



Y,K>Y,hi (113) 



Hence the supremum cannot be taken on the interior. We conclude that at least 
some hf_^_i + hf have to be one. There is then a minimal i for which this holds. 
We can convince ourselves that the process can now be repeated for the substring 
ending at i + l. Repeating the argument we can fill up all the entries up to place 
i + 1 with the condition hf^-^ + hf = 1 and proceeding now upwards we end up 
with 
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Lemma 4.19 The above supremum is assumed within the subset 

hl<l,hl + hl = l,... , hl_, + hl = l,hl<l (114) 

This concludes the first step. 

In the second step we calculate sup |/(0) — f{n) \ on this restricted set. From 
the above we now have the constraint: 

hl<l,hl = l-hl,hl = hl,hl = l-hl,... ,hl = l- h\ or hi (115) 

depending on n being even or uneven. This yields 

1 for n = 1 



sup 1/(0) - /(n)| = <j (n/2) ■ sup(/ii + - hi) = {n/2) ■ ^2 for n even 

^sup([n/2] ■ {hi + \/l — hi) + hi) for n uneven 

(116) 



In the even case the rhs can be written as v/nV2 = v/K/2]. In the uneven case 
we get by differentiating the rhs and setting it to zero: 



/zr^ = Aj^i + , vi - {hrn' = v vi + (n?) 

with An = 1 + l/[n/2]. We see that for increasing n both terms approach 1 / a/2, 
the result in the even case. Furthermore we see that the distance is monotonely 
increasing with n as should be the case for a distance. This yields in the uneven 
case 

, , (\n/2] + l)An+\n/2] , , 

distc{<d, n) = '±±-L^=^^±A (118) 

which is a little bit nasty. Both expressions can however be written in a more 
elegant and unified way (this was a conjecture by W.Kunhardt, inferred from 
numerical examples). For n uneven a short calculation yields 

[n^/2] = (n^ - l)/2 = 1/2 ■ (n - l){n + 1) = 2[n/2]([n/2] + 1) (119) 

(with the floor- , ceiling-notaXion the expressions would become even more ele- 



gant). With the help of the latter formula the rhs in ( |118|) can be transformed 
into 

i[n/2] + l)An + [n/2] 



VhV2] + 1 (120) 
Conclusion 4.20 For the one- dimensional undirected lattice we have 



N 1 \fWl% for n even 
dzstc{0,n) = {^ \ / ^ (121) 

y [n^/2J + 1 for n uneven 
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Remark 4.21 Again comparing the exact result with our lower bound, we find 
for n even: 



distciO, n) > (2)^^/2 . ^ ^ ^2)^/2 . (122) 
that is, the lower bound is again saturated. For n uneven we have instead: 

distciO, n) = (2)-i/2 . (^2 ^ > ^2)-^/2 . ^ ^^23) 

Example 3: The directed lattice 

The vertices in are denoted by (i, j) or {x,y). The edges point from {i,j) to 
{i + 1, j) and {i,j + 1); hence, ff"* = 2. As the system is translation invariant, 
it suffices to calculate the Connes distance between nodes (0, 0) and (x, y) with 
x,y > 0. 

For nodes lying on the same parallel to the x-, y-axis, respectively, we have 

distc{n,n) = d{n,n') (124) 

For X or y = 0, there is only one minimal path, connecting (0,0) and {x,y). 
Therefore, lemma |4.8| does not apply. For, say, y = 0, we choose the following 
admissible function: 

f{x,y):=x for all y (125) 

We have 

\f{x, 0) - /(O, 0) I = |x| = d{{0, 0), (x, 0)) (126) 

and can conclude 

distciiO, 0), (x, 0)) = d{{0, 0), (x, 0)) (127) 

on Z^. The same holds for the y-axis. 

For nodes with both x,y 0, we have more than one minimal path connecting 
(0,0) and {x,y). Our lemma then shows that, necessarily, 

distc{n,n') < d{n,n') (128) 

More detailed estimates will be given below. 

If we try to really calculate the Connes distance on for points in general 
position, the optimization problem becomes quite involved and we will only pro- 
vide some estimates. The reason is that the constraint equations are of a quite 
non-local nature (compared to the simpler undirected Z^-lattice) and that in gen- 
eral several minimal paths do exist which make the continuation problem quite 
intricate. 
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It is easy to see that the canonical graph distance between the points (0, 0) and 
(x,y) is |x| + \y\ and that all minimal paths have the same length. With x, ?/ > 0, 
we conjecture (without giving a proof) that it suffices to restrict the variation 
to admissible functions with positive jumps in the positive x— , |/— directions and 
that we can set /(O, 0) = 0. A particular minimal path consists of x steps in the 
x-direction followed by y steps in the y-direction. We denote (for convenience) 
the jumps along the a;—,?/— axis, respectively, by 

hio := /(^,0) - /(z - 1,0) > /lo, := /(0,j) - /(0,j - 1) > (129) 

The optimization problem now reads: 

Problem 4.22 Find sup ^X]i=i ^io + X]J=i ^oj j under the constraints imposed 
by the admissibility of the corresponding function, f . Note however that the 
constraints must hold on the full lattice. 

From our general result in theorem [4.15| we know that 

dzstcm 0), {x, y)) > (2)-i/2 . (3; + ^) (130) 

We can construct an admissible function, /, which fulfills 

1/(0,0) -/(x,i/)| = (2)-i/2-(x + y) (131) 

This can be achieved by setting /(0,0) = and by choosing all x—,?/— jumps 
equal to a with 2a^ = 1 ^ a = (2)~^/^. This yields the above result. 

The question is, whether this is already the supremum over the set of admis- 
sible functions. We will show, that this is not the case by providing an other 
admissible function yielding a bigger value. We choose an admissible function 
with x-jumps equal to a and y-jumps equal to b with + 6^ = 1. The admissible 
function reads 

f{x,y)=ax + by (132) 
The function / takes a stationary value at 

a = ix'/ix' + y')Y'' b = iy'/ix' + y')fl' (133) 

yielding the value 

f{x,y) = {x' + y'Y/^ (134) 

Assuming, for example, that x 7^ y, it follows that 

/(x, yf - ((2)-V2 .{x + y)f= (1/2) . [x^ + y^) - x ■ y = (1/2) . {x - yf > Q 

(135) 

In other words, we see 

Observation 4.23 For the directed W?-lattice and x ^ y (x,y > 0) we have the 
estimate 

d{{0, 0), (x, y)) > dtstciiO, 0), (x, y)) > (x^ + yY' > (2)"'/'(x + y) (136) 
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